%% 本算例：(########算例说明文件头###########)

% lamda=c*T=2*pi/w
% 例如：lamda=200, w=0.01*pi

% 本算例：无初始条件*，只有边界条件（位移单峰波）**
% 对应abaqus中，没有distribution，只有amplitude

%% 初始化：建立work
clear
clc 

workname=input('Input the job name:','s');

dirpath=strcat('dataSaved/',workname,'/');
mkdir(dirpath);
workname_path=strcat(dirpath,workname);
work_log=strcat(workname_path,'.log');

fid=fopen(work_log,'w');
fprintf(fid,'[ODCpre] Work Name: %s',workname);

time_start=datestr(now);
fprintf(fid,'\n[ODCpre] Start: %s\n\n',time_start);

YES=1;
NO=0;

%% 1. 设置：模型参数

ctn=1;  % 密集指标（越大越连续体，1表示标准离散化，越小越离散）
L=1000; % 模拟连续杆件的长度

cm=0; % 瑞利阻尼a0
ck=0; % 瑞利阻尼a1

fprintf(fid,'[ODCpre] ctn = %d\n',ctn);
fprintf(fid,'[ODCpre] L = %d\n',L);
fprintf(fid,'[ODCpre] Damping %0.3f*M+%0.3f*K\n\n',cm,ck);

%% 2. 设置：分析步

dt=0.001;	%物理时间间隔
Tt=500;     %物理时间总时长

outdt=1; % 场变量输出时间间隔 
dframe=outdt/dt;	% 场变量输出频率，每隔outdt时间输出一次，也就是每隔dframeA帧输出一次

fprintf(fid,'[ODCpre] dt = %f\n',dt);
fprintf(fid,'[ODCpre] Tt = %0.2f\n',Tt);
fprintf(fid,'[ODCpre] outdt = %0.2f\n',outdt);

VV=101;
Newmark=102;
option=101;

fprintf(fid,'\n[ODCpre] option = %d\n',option);

%% 3. 设置：初始条件*

%   SetIniCondition.m函数内设定
fprintf(fid,'\n[ODCpre] Not initial conditions yet.\n');

%% 4. 设置：边界条件**

w=0.01*pi;
t_cut=100;

fprintf(fid,'[ODCpre] w = %f\n',w);
fprintf(fid,'[ODCpre] t_cut = %d\n',t_cut);
    
 %% 其他设置
 
switchINFO=YES;	%信息输出设置，默认输出

fprintf(fid,'\n[ODCpre] PreProcess over.\n');
fprintf(fid,'\n=================================================\n');
fprintf('[ODCpre] PreProcess over.\n');